Adiabatic Approximation from a Renormalization Group Inspired Method 
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Two related renormalization group schemes for Hamiltonians have been independently proposed in the last years 
by Wegner |Q and Glazek and Wilson Q in the context of condensed matter and field theory, respectively. These 
nonperturbative methods are useful for diagonalizing or renormalizing a given Hamiltonian and are based on a 
continuous sequence of (infinitesimal) unitary transformations applied to the system Hamiltonian. The key point 
in these approaches is that the unitary transformations are chosen in such a way that off-diagonal matrix elements 
becomes smaller at each step of the sequence. The outset of this procedure is a set of flow equations whose solution 
yields (in the best case) a (almost) diagonalized Hamiltonian. These methods have been succesfully applied to various 
problems in condensed matter and light front field theory J3|. 

In this letter we take a look at the origins of this kind of renormalization transformations, i.e., the work of Frohlich Q] 
(see also ||). He gave an explanation for the effective interaction between electrons in BCS theory of superconductivity 
by eliminating electron-phonon interactions using a unitary transformation, which can be viewed as the "one-step" 
formulation of the method mentioned above. We show how Frohlich method can be used as an approach to the usual 
adiabatic approximation used in quantum optics. To that end we will use as a toy model the Dicke Hamiltonian Q| 
which describes a collection of two-level atoms interacting with a quantum field. 

In the method of Frohlich a similarity transformation is performed on the Hamiltonian such that the transformed 
Hamiltonian has the (troublesome) off-diagonal elements equal to zero. Since similarity transformations preserve 
eigenvalues, the transformed Hamiltonian has the same spectrum as the original one. Thus the new Hamiltonian is 
given by 

H v = UHUl (1) 

Let us write the unitary operator as 

U = e~ m (2) 

where SI is defined as a perturbative series in a coupling constant. If we separate the Hamiltonian into its free and 
interacting parts we have 

Hu = H + {Hi + i[H , SI]) + i[Hi,(l] - -[{H ,Si\,Si\ + ■■■ (3) 

We must look for an operator SI such that it eliminates the effect of the Hamiltonian Hj to lowest order, that 
is, such that the transformed Hamiltionian is completely diagonal to that order. This means, we must impose the 
condition 

Hi + i[H ,Sl] = (4) 
which in turn implies that when this condition is satisfied, Hjj can be written as 

H u = H Q + ~[H I ,Sl} + --- (5) 

To show how the adiabatic approximation is obtained from it we use as a toy model the Dicke Hamiltonian 

H = LUoa^a + u>iS z + g (aS+ + a i S-) (6) 

where a and a' are the usual field operators and the collective atomic operators satisfy 

[S+,S-] = 2S Z , [S„S±]=±S± (7) 

In the limit of large detuning, uj\ luq, we can adiabatically eliminate the transitions among different eigenstates 
of S z . That is, we want to decouple the low frequency modes from the high frequency modes (this can be done at 
least perturbatively by diagonalizing with a suitable unitary operator). By replacing Ho and Hi into eqs. (||) and (^) 
we have 

[uoa ! a + iviS z , SI] — ig {aS+ + a i S-} (8) 

From the structure of this equation and the algebra of the operators involved, one can easily infer that O must be 
necessarily of the form 

Sl = g{aaS++ Pa^S-} +SI' (9) 
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where CI' is any operator that commutes with the free Hamiltonian Hq. That is, the solution to eq. (Q) is always of 
the form 

ci = ci + ci' (10) 

The fact that CI' remains undetermined is harmless. In effect, we have 

H v = e -x*-in H ja'+i(i 

but we can always perform a unitary transformation on Hjj such thatQ 

H' v = e- m e- in 'He in 'e m . (12) 
which means that an unknown unitary transfomation is acting on the original Hamiltonian H: 

H' = e.- iU 'He iU ' (13) 

Since a unitary transformation does not change the eigenvalues of the Hamiltonian, we can always make the choice of 
working either with H' or H. By choosing to work with H is equivalent in practice as imposing the condition CI' = 0. 
We do this in the following. 

Therefore, after replacing eq. (ph into eq. (@) we get 



Cl = i^- {aS+ - a) S-} (14) 



where A — lu\ — luq. 

Also, the 0(g 2 ) contribution 



is readily obtained from ffl) and (B 



H^ = -[H!,Cl} (15) 
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#(2) = y {2S z a^a + S+S- } (16) 
which implies that the correct total Hamiltonian to this order in perturbation theory is 

2 2 

H LI = woo+o + lo!S z + 2^aS z + ^ (S 2 - S 2 Z + S z ) (17) 

This Hamiltonian can also be derived by using density matrix techniques in the adiabatic approximation but the 
calculation is too much involved Is]. It has also been derived by using unitary transformations in H, but the way of 
obtaining it was based on the calculation done in and then performing a guess. In that way they obtained 

2 2 2 

H = Looa^a + uj.S, + 2^aS z + |- (S 2 - S 2 + S z ) + ^ (S 2 + + Si) (18) 

where the last term had to be discarded by making a further approximation. 

One could have also used a sequence of discrete unitary transformations, as done e.g. in H. After a lengthy 
calculation one gets 

H v = uo^a + lj^z + (^ + ^ + ^1) {2S z a)a + S+S-} + ^ {aS+ + a^S\ (19) 

UJl Luf LU L { Ujf 1 ' 

where the last term must be removed by another unitary transformation and we recognize in the remaining 0(g 2 ) 
term part of the geometric series corresponding to the expansion of <? 2 /A. 



1 This can be easily proved by following the standard steps one makes when deriving the Baker-Campbell-HausdorrT or related 
(e.g Zassenhaus) formulas. 
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In this letter we have shown that it is possible to use a renormalization group inspired method to derive Hamiltonians 
in the adiabatic, large detuning, approximation in a straightforward and systematic way, avoiding the appearance of 
extra terms that should not be there. It must be stressed that although the model we used to show the method is 
simple, the procedure can be implemented with no difficulties in more realistic models. The crucial step will always 
be to find the solution to the operator equations that appear in the procedure, like eq. (^) above. This is done by 
inspection, which in practice it is not difficult to do. 
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